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We propose a scheme for extracting an effective three-body interaction originating from a two- 
nucleon interaction. This is based on the Q-box method of Kuo and coUaborators, where folded 
diagrams are obtained by differentiating a sum of non-folded diagrams with respect to the starting 
energy. To gain insight we have studied several examples using the Lipkin model where the per- 
turbative approach can be compared with exact results. Numerically the three-body interactions 
can be significant and in a matrix example good accuracy was not obtained simultaneously for both 
eigenvalues with two-body interactions alone. 

PACS numbers: 21.60.Cs, 21.30.Fe, 24.10.Cn, 21.60.Fw 



I. INTRODUCTION 

The role of effective three-body forces in shell-model 
and nuclear structure calculations of medium-heavy and 
heavy nuclei is still an unsettled question. The question 
has received scattered attention over the years, but there 
has been no decisive studyyet yielding a unique answer, 
see for example Refs. [1 H H 0. Most of the studies 
seem to indicate that effective three-body forces play a 
rather modest role for nuclear spectra, although there 
is a possibility that binding energies may be more sub- 
stantially affected. As pointed out by Zuker Isj, one is 
not able to obtain simultaneously a good reproduction of 
both the excitation spectra and the binding energy with 
effective interactions based on two-body nucleon-nucleon 
interactions only, unless one fits an effective interaction 
to reproduce selected data, see for example the recent 
work of Otsuka et al. in Refs. 0,0- To give an exam- 
ple, effective interactions derived from two-body nucleon- 
nucleon interactions which fit nucleon-nucleon scattering 
data, are not able to reproduce the well-known shell clo- 
sure in '^^Ca nor the excitation spectra of ^^Ca and "^^Ca 
0, 0, 0] ■ For the chain of oxygen isotopes one may even 
reach the conclusion that ^^O is a bound nuclear system, 
unless the interaction is fitted to reproduce selected prop- 
erties such as binding energies and spectra of known and 
stable nuclei in the \sOd shell @. 

On the other hand, theoretical interactions which have 
not been fitted do rather well in reproducing nuclear 
properties other than binding energies and shell clo- 
sures. Examples are nuclear s pec tra around A ^ 100 
and A - 132 [IS ini mill HI [13, which come out sur- 
prisingly well in view of the complexity of the many-body 
problem, ft has therefore been speculated that some of 
the abovementioned deficiencies could be ascribed to the 
lack of three-body interactions, as seen in recent Monte 
Carlo and no-core shell- model calculations of light nuclei 

with A < 16 ElliLmia, 201). 

In light nuclei the role of three-body forces seems 



rather well established. To reproduce the binding energy 
of the triton, for example, it is necessary to include a 
three-nucleon interaction in addition to the two-nucleon 
interaction. The situation in medium-heavy and heavy 
nuclei is however unclear, although similar arguments, 
based on an analysis of Op, IsOd and IpO/ nuclei by Zuker 
13 lend support to the need for three-body interactions 
in heavier nuclei as well. 

In contrast to the case for light nuclei, none of the 
present calculations of effective interactions for medium 
and heavy nuclei are sufficiently complete, however, to 
enable one to draw definite conclusions. For nuclei with 
A > 16, one relies on shell-model analyses combined with 
appropriate effective interactions for a selected model 
space. Such interactions are typically based on perturba- 
tive many-body methods and because high-order calcu- 
lations of the effective three-body contribution are pro- 
hibitively complicated, one is forced to limit oneself to 
low-order calculations. There are however classes of di- 
agrams which can be summed to infinite order. One of 
these classes is the set of folded diagrams. The latter 
arise due to the removal of the dependence on the ex- 
act model-space energy in the Brillouin-Wigner pertur- 
bation expansion. Through the Q-box formulation and 
its derivatives 0, |^ , this set of diagrams can easily be 
summed up. With few exceptions Q], three-body folded 
diagrams have not been included. There are indications, 
however, that folded three-body diagrams may be impor- 
tant, even if non- folded ones are small. 

There is obviously a need to obtain a clearer picture of 
the role of effective three-body forces in heavier nuclei. 
Part of our motivation for studying the contribution of 
effective three-body forces comes from large-scale shell- 
model calculation of the entire range of Sn isotopes from 
A — 100 to A = 132, where the spectra are well repro- 
duced whereas the binding energies are strongly overes- 
timated. In fact, it was found that a small repulsive 
monopole contribution to the effective two-body interac- 
tion could cure this problem. The origin of such a term 
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is unknown, but it is speculated that it might be simulat- 
ing an effective three-body interaction. Indeed Talmi 
has shown that the binding energies of the entire range of 
Sn isotopes could be well fitted by an effective two-body 
plus an effective three-body interaction. In spite of the 
fact that both the Green's function Monte Carlo and no- 
core shell-model approaches offer benchmark calculations 
of light nuclei with A<16ll6|,ll7|,ll8|,lll|23, and that 
coupled cluster approaches can presently extend the re- 
gion of ab initio calculations to A ^ 40 — 56 , features 
of heavier nuclei will most likely be studied within the 
framework of the shell model and appropriately defined 
effective interactions. Thus we need to have a controlled 
approach which allows us to discern contributions from 
various many-body terms, such as effective three-body 
forces. Effective three-body forces may originate from 
three-nucleon forces as well as from two-nucleon forces. 
In order to be able to draw definite conclusions, we be- 
lieve that it is necessary to study one effect at a time. 
Thus, before introducing a three-nucleon interaction, we 
would study the contribution from effective three-body 
forces originating from the two-nucleon interaction. This 
is the topic of the present work. 

The paper is organized as follows. In Sec. |n]we pro- 
pose an algorithm for evaluating the effective three-body 
interaction in the convenient Q-box formalism introduced 
by Kuo and collaborators We prepare ourselves 

for this procedure by briefly summarizing the method for 
obtaining the effective two-body interaction. Because the 
method can be rather involved in realistic cases, we apply 
it in Sec. Illll to a simple model, namely the Lipkin model 
|23 |. for which exact results can easily be obtained, see 
for example Ref. j2^. This is instructive both as regards 
the formalism and the numerical results. In the final Sec. 

we draw conclusions and point to the need for more 
realistic calculations. 



II. THE EFFECTIVE ONE-, TWO- AND 
THREE-BODY INTERACTIONS 

The Hamiltonian acting in the complete Hilbert space 
(usually infinite) consists of an unperturbed one-body 
part, Hq, and a perturbation V, namely H ~ Hq + V. 
The goal is to obtain an effective interaction, V, such 
that H = Hq + V acting in a chosen model or valence 
space yields a set of eigenvalues which are identical to (a 
subset of) those of the complete problem. To this end 
Kuo and collaborators introduced the Q-box, see Refs. 
0, which, in principle, contains all non- folded di- 
agrams which are attached to the valence particle lines. 
Unlinked pieces which refer to excitations of a closed-shell 
core are absent since the energy is calculated relative to 
the closed-shell ground state [2g. The effective interac- 
tion also requires that folded diagrams be included to all 
orders and this can be formally written 



where 



V{n) = Qiio) + J2 



1 d" 



■{V{n-l)y 



(2) 



with 1^(0) = Q and matrix multiplication implied. Here 
uj is the unperturbed energy of the initial state. In or- 
der to perform calculations with N particles in the model 
space it is necessary to separate the effective interaction 
V into its one-, two- and three-body components (in prin- 
ciple up to A'^-body interactions are needed, but we shall 
not go beyond three-body) . This brief sketch establishes 
our notation; a complete discussion is to be found in Refs. 
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V = lim V{n) 



(1) 



(a) (b) (c) 

FIG. 1: Two-body diagrams discussed in the text. 

Consider first the case where our model space consists 
of a single particle, N = 1; we assume that a given total 
angular momentum j is represented only once. Then Eq. 

can be solved for each state separately, yielding the 
complete folded series for the one-body effective interac- 
tion V = The one-body effective Hamiltonian is 
thus Ho + V^^\ 

Now consider the case where we have two particles in 
the model space, N = 2. The Q-box again consists of 
all non- folded diagrams. It will include diagrams such 
that the interactions link the two valence lines together 
as indicated schematically in Fig. Q^a), diagrams where 
one line is non-interacting as in Fig. Q^b) and also di- 
agrams where the interactions on each line are separate 
so that the diagram is valence unlinked as in Fig. ^c). 
All of these components need to be included so as to 
generate the cross terms in the full folded effective inter- 
action of Eqs. ^ and ((2Jl. Valence unlinked diagrams 
are removed by the folding process so that the effective 
interaction is completely linked [2^ . Now because of the 
one- plus two-body nature of the Q-box, the result for V 
from Eq. ^ will contain linked, folded two-body com- 
ponents which are schematically of the form of Fig. QJa) 
plus linked, folded one-body components where one line 
is non-interacting, as in Fig. H^b). The one-body piece 
V^^'^ is already determined, as indicated above, so the 
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FIG. 2: Three-body diagrams discussed in the text. 

purely two-body interaction y^^^ can be obtained from 

{de\V\ab)2 = {de\V^^^\ab)+5daSeb{(a\V^'^\a) + (b\V'^'^\b)y 

(3) 

The subscript on the left indicates the value of N for 



which Eq. is solved. We assume here that the ini- 
tial and final states in Fig. ^a) are in a standard order 
so that e = a implies d ^ b. Here the two-body ma- 
trix elements are "antisymmetrized" , i.e {de\V\ab) im- 
plies {de\V\ab — ba). They are not coupled to a total 
angular momentum or isospin since our shell model cal- 
culations are carried out in the m-scheme. 

Turning to the more interesting three-body case (N — 
3), consider the Q-box. It contains diagrams such that 
the interactions link the three valence lines together as 
indicated schematically in Fig. EJa), diagrams where two 
of the lines are linked (Fig. El^b)) or where only one line 
is interacting as in Fig. [2Ic). The valence unlinked dia- 
grams can involve a two-body and a one-body piece (Fig. 
|2Id)) or two or three one-body pieces (Figs, ^e) and (f)). 
Notice that there are no valence unlinked diagrams if, as 
is often the case, one-body interactions are excluded. In 
general all of the components illustrated in Fig. |21need to 
be included in generating the full folded, valence-linked 
effective interaction of Eqs. ^ and Q. The interaction 
has a three-body contribution where the three lines are 
linked together, schematically in the form of Fig. El^a), 
two-body interactions acting on all possible pairs of lines 
(Fig. |2Ib)) and one-body interactions acting each line 
separately (Fig. [2Ic)). Using Eq. (PJ we can disentangle 
the desired V"''^^ from 



{def\V\abc)3 = {def\V^-^^abc) + 5 fc{de\V\ab)2 - 5fb{de\V\ac)2 + S fa{de\V\bc)2 + 5eb{df\V\ac)2 
-Sec{df\V\ab)2 - S,a{df\V\bc)2 + Sda{ef\V\bc)2 - Sdb{ef\V\ac)2 + Sdc{ef\V\ab)2 

-SdaSeb6fc{{a\V^'^\a) + {b\V^'^\b) + {c\V^'^\c)) . (4) 



Again the subscript on the left indicates the value of N 
for which Eq. ^ is solved. Of course for a given matrix 
element only a few of the delta functions will be nonzero. 
Since the inital and final states are in a chosen standard 
order the one-body term only contributes when d = a, 
e = b and f = c. Here the uncoupled three-body matrix 
elements are "antisymmetrized", i.e {def\V\abc) imphes 

{def\V\abc — acb + cab — bac + bca — cba) , 



labelling in order for particle numbers 1, 2 and 3. Equa- 
tion 10} has one-body components in {V)2 as well as a 
compensating explicit term. By using Eq. Q the equa- 
tion can be cast just in terms of one-, two-, and three- 
body effective interactions: 



{def\V\abc)3 = (de/lF^^VM + ^/c(rfe|F(2) |a&> - '^/bidel^^^Vc) + '5/a(de|l^(2) |^ ^ 5^^^^j-|y(2) |^^^ 
-d,M\V^^^\(ib)-S,a{df\V^^Hbc)+dda{ef\V^'^bc)^ddb{ef\V^^^\ac)+Sd,{en^^^^ 
+6da6eb5fc{{a\V^'^\a) + {b\V^'^ \b) + {c\v'-'^c)) . (5) 

I 



This just spells out explicitly the fact that the complete effective interaction requires all possible valence-linked 
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contributions. 

Thus, given the Q-box in the one, two and three parti- 
cle systems we can obtain the one-, two- and three-body 
effective interactions. Obviously it would be possible to 
generalize this procedure to obtain four-body and higher 
interactions. However, since almost all calculations have 
stopped at the two-body level, it is sensible to investigate 
just three-body interactions. Here we attempt to get a 
feel for their impact in a simple model situation. 



III. APPLICATION TO THE TWO-LEVEL 
LIPKIN MODEL 



The Lipkin model |2J| consists of two single-particle 
levels labelled by ct = — and each of which has a 
degeneracy p. We write the Hamiltonian 



J 



H 

Ho 



= Ha + V 

= KE- 

pa 



where 



and 



pp'a 



■ a P~ 



t t 

pa^p'a'^P' 



pp' a 



t t 
^po^p' - 



■ a P 



,(6) 



Here Hq is the unperturbed Hamiltonian with single par- 
ticle energies ±5$. The two-body interaction, V, has 
three terms. The interaction V acts between a pair of 
particles with parallel spins and changes the spins from 
-|--|- to , or vice versa. The interaction is a spin- 
exchange interaction and U, which was not present in the 
original model flips the spin of one particle. It is of 
interest to note that the interaction does not change the 
value of the degeneracy labels pp' . 

Since each particle has only two possible states, the 
use of the quasi-spin formulation was suggested by Lip- 
kin et al. [2j|. The quasi-spin operators obey angular 
momentum commutation relations and are defined by 

Jz — 2 ^ ^ '^Qpo-Qpg- J J+ = ^ ] ^p+^p- I ^ ^ ] '^p~'^p+ ■ 

per p p 

(7) 

The Hamiltonian can then be compactly expressed in the 
form 

H = eJ^ + il^(J+ + J') + iVF(J+J_ + ,/_,/+ -n) 
+ i{/(J+ + J_)(n-1) , (8) 



where the number operator n = X^po- '^po-^po-- The op- 
erator = i(J+,/_ + J-J+) + Jl commutes with the 
Hamiltonian so the Hamiltonian matrix breaks up into 
submatrices of dimension 2 J -I- 1, each associated with 
different values of J; for a given number of particles N 
the largest angular momentum corresponds to J = \N . 

Here we need the case iV = 3 and for our purposes it 
is sufficient to consider J — ^. We denote the basis state 
for three particles by 



|pi±,p2±,P3±) = api±"p2±43±|0> ' (9) 



where pi refers to the degeneracy label and |0) is the 
vacuum state. Then the basis states for the = | matrix 
are: 





\J = 


2'Jz — 


-1) 


= \P1~ 


,P2-,P3-) , 










\J = 




-h) 


= {\P1 


-,P2-,J33+> 


+ \Pl-,P2 + ,P3-) 


+ \Pl~ 


\^:P2-,P3-)}/V3 , 




\J^ 


2'Jz — 


+ \) 


= {\P1 


~,P2 + ,P3+) 


+ \Pl + ,P2-,P3 + ) 


+ \Pl~ 


^,P2 + ,P3-)}/V3 , 


|4)^ 


\J = 


\Jz = 


+ 1) 


= bi+ 


,P2 + ,P3 + ) ■ 









In this basis the Hamiltonian matrix is 



The Case f/ = 



V 



2U 



2U 



\ 

k J 



(11) We consider first the simplest case where U ~ so that 

the matrix Hll|l splits into two 2x2 matrices. Consider 
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(a) 
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(c) 



FIG. 3: Third-order (a) and fourth-order ((b) and (c)) three- 
body diagrams discussed in the text. Diagram (b) is a folded 
diagram drawn in unfolded form for clarity. 



the matrix formed by states |1) and |3): 



\-/3V in 

whose eigenvalues are 



2W 



1± 



(12) 



(13) 



In principle both eigenvalues are given by an exact treat- 
ment of the perturbation series 28] since Eq. ^ is sim- 
ply a rearrangement of the exact eigenvalue equation. In 
practice approximations imply that only one eigenvalue 
is obtained and we will focus on A_ using the model 
space state |1). 



Now directly from the matrix H12|) . or by explicitly 
drawing diagrams, the second order Q-box is 3V'^ / (—2^). 
Clearly the W interaction can be incorporated to all or- 
ders yielding an exact Q-box of 3V^/{-2^ - 2W). The 
folded diagrams can then be included by differentiation 
as indicated in Eq. |(5J). In the approximation that we 
stop at order the interaction is 



(^)3 = - 



3V' 



^ MO ■ (14) 

In order to disentangle the three-body interaction we 
need the two-body interaction for the state \pi — ,p2—). 
Since this is only coupled to \pi+,p2+) by the Hamil- 

V ' 
V ^ 

use we note that the exact lowest eigenvalue is A!l^ = 
-^y^l + -p- Now the exact Q-box is V'^/{^ 
rying out the folding to order as before 



tonian, the necessary matrix is 



For later 

(2) 



T/2 y4 yG 

(^h = -^ + ^-^^/.(0. (15) 



The one particle wavefunction is ap_ |0) with energy — 
given by the unperturbed Hamiltonian Hq, so that the 
one-body interaction V'^^^ is zero. Thus, using Eq. (0)) 
or , and noting that because of the degeneracy label p 
there will only be three contributions from the two-body 
terms 



~2e~ 



4^3 



2W^) 



W^W^ + 2W'^) . (16) 



r 



Here we have expanded to sixth order in the interaction. 
Notice that there is no second-order contribution. This 
result can also be obtained by evaluating diagrams, for 
brevity we do not go beyond fourth order. In third order 
the diagram of Fig. Ola) is required, while in fourth order 
the folded diagram of Fig. O^b), drawn in unfolded form, 
and the non-folded diagram (c) are needed; by counting 
the number of independent diagrams of each type the 
numerical factors in Eq. (|16() are obtained. If W is set to 
zero in Eq. H16(l the remaining three-body contributions 
come from folded diagrams. Fig. Otb) and higher orders. 
The coefficients of the order V'^ and terms in Eq. 
(|16|l are larger than three times the corresponding values 
in Eq. H15|l . This supports the suggestion that even 
if non-folded three-body diagrams are small, three-body 
contributions arising from folding could be important. 



We can, of course, obtain the exact three-body inter- 
action for this simple Hamiltonian using the exact eigen- 
values given. As we have remarked, the one-body inter- 
action yd) is zero. Then, subtracting the unperturbed 
energy, the two-body effective interaction is 




(17) 



The desired exact three-body interaction is therefore 



(W3)) = AL'^-3(y(2)) + 3^ 



(18) 



where the last term removes the contribution from _ffo in 



A 



(3) 



Expanding this expression to sixth order of course 



yields Eq. (|16f) again. Notice that since V^'^^ depends 
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(a) 



(b) 



(c) 



FIG. 4: Fourth-order three-body diagrams discussed in the 
text. Diagram (c) is a folded diagram drawn in unfolded form 
for clarity and the dot represents a one-body insertion. 



on W while y*^^-' does not, suitable adjustment of this 
parameter could yield a three-body interaction of mag- 
nitude comparable to the two-body interaction. For in- 
stance, in the limit Vl^^ ^ cxd we have ^ _3\/(2) 
since the eigenvalue is simply — |^ in this case. 

It is instructive to make the model a little more elabo- 
rate by changing the parameter ^ of the one-body Hamil- 
tonian to ^-f e. The first term is included in Hq as before 
and e is treated as a perturbation. The matrix (|12|l then 
becomes 



From the matrix directly, or by drawing diagrams, the 
Q-box through second order contains both one- and two- 
body parts, namely — |e — 3F^/(2^). Starting in third 
order there will be diagrams with e insertions. Since the 
intermediate state always involves two -t- particles and 
one — particle, the net insertion is +ie (note that va- 
lence unlinked contributions of the form of Fig. |2Id) 
occur here). If these insertions, together with the W in- 
teractions are summed to all orders one obtains — |e — 
3V^/{2^ + 2W+^e), a result that can be obtained imme- 
diately from the matrix (|19|l . Now an e insertion on each 
of the three — particles of the model space state |1) gives 
a net contribution of — |e. If these insertions are folded 
out to aU orders one obtains -|e - 3^ V(2C + 2W + 2e), 
in the process removing valence-unlinked diagrams. This 
result again follows immediately from the matrix H19|l . 
If this modified Q-hox is used in Eq. Q to generate 
all the additional folded diagrams required we obtain 
{y)3 — — |e + fsiS, + e)- A similar procedure in the two- 
body case gives {V)2 — — e + /2(C + f)- the one-body 
case (1/(1)) = -if. Then using Eq. ^ 



-3£_ 3 



2W 



(19) 



(F(3)) = (1^)3-3(^)2 + 3(^1)) 

3T/2M/ 3T/2 oT/2 

= — — -1- —:r(y^ - 2M^^ - 4:eW) + -—(4W^ - 9V^W - 6V^e + 12W^e + l2We^) 
2^2 4^3 ^ ^ 8^4 ^ ' 

-r{2V'^ - W^W^ + 2W^ - mV^We + 8W^e - 6V^e^ + UW^e^ + 8We^) , 



(20) 



to sixth order in the interaction. Obviously the terms 
which do not involve e are the same as in Eq. No- 
tice that terms of the form T/2giy^n+i necessarily belong 
to the two-body interaction and do not occur here. The 
fourth order term involving e is generated from the dia- 
grams of Fig. Efa) and (b), as well as the folded diagram 
(c). Here V'-^^ is represented by a dot. Each diagram 
stands for several different diagrams of the same general 
structure and the coefficient in Eq. H2(J|I is obtained by 
accounting for all of them. The result in Eq. (|20|) can 
also be obtained by a power series expansion of the exact 
result which is the analogue of Eq. (|18|l . namely 

(1/(3)) = aL'^ - 3(t/(2)) - 3(^(1') + |e , (21) 

with the replacement ^ ^ ^ + e being made in the first 
two terms on the right. 



B. The Case [7/0 



Here we take the model space to consist of states |1) 
and 1 2) of Eq. H1U|I and require our effective interaction 
to reproduce two of the eigenvalues of the matrix Hll|l . 
It is simplest to work with a degenerate model space and 
we choose the unperturbed energy to be — The re- 
mainder of the diagonal (2, 2) contribution in (|ll|l is then 
included as a perturbation. The lowest order Q-box of 
a given type is easily written down as discussed in Sec. 
nil Al Here we also need to account for multiple scatter- 
ings back and forth between states |3) and |4) which is 
easily done by summing the infinite series. Thus we can 
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> 
c 

9) 




1 2 3 4 5 
fold fold fold fold fold fold exact 



FIG. 5: Comparison of the eigenvalues (in units of ^) of the 
effective interaction with n folds to the exact result. For 
the solid lines two- and three-body interactions are included, 
while for the dashed lines only two-body interactions are con- 
sidered. 



obtain the exact Q-box: 



(lIQIl) 



(2|Q|1) 



-2£_-2W + U^/^ 
(l|g|2) = V3U 



2V3UV 



(2|g|2) =e + 2VK 



2e(e + w)-u^ 



-2^-2W + U^/^ 
WU^ 

-^ + Uy{2^ + 2W)^ 2^{^ + W)-m ' ^ ' 

Insertion of this in Eq. ^ will yield an effective in- 
teraction which reproduces two of the exact eigenvalues. 
This is illustrated by the solid lines in Fig. [S] We work 
in units of ^ and have chosen the dimensionless ratios 
V/^ = VK/C = -0.4 and U/^ ^ -0.15. The label fold 
in Fig. |S1 implies that the solid lines are obtained from 
the Q-box alone, while the labelling n fold implies that 
the summation over m in Eq. (j^J extends to m = n. As 
n is increased so that more folded diagrams are included 
the eigenvalues of the effective interaction converge to the 
exact result. The deviations are less than 0.2% by the 
time n = 5. 

This analysis necessarily includes effective three-body 
interactions and we would like to assess the effect of their 
removal. To that end consider first the two-body prob- 
lem. We only need to consider the J = 1 matrix which 
is 



U/V2 V 
U/V2 W U/V2 
V U/V2 C 



(23) 




FIG. 6: A three-body diagram discussed in the text. 



in the basis 



->}/V2, 
(24) 

Taking the model space to comprise states \a) and \b) 
and choosing degenerate unperturbed energies of — ^, the 
exact Q-box is 



a) ^ 


\J = 


IJz = -1) 


= \pl- 


-,P2-) , 


\b) ^ 


\J = 


IJz = 0) = 


{IPI- 




\c) ^ 


\J = 


IJz = +1) 


= \PI^ 





(«IQI«> = 



2e ' 



m\a) 

{b\Q\h) 



{a\Q\b) 



v_ _ vu 

' V2 2V2( ' 

El 

4e ■ 



(25) 



When used in Eq. ^ this will give an effective interac- 
tion which yields two of the exact two-body eigenvalues. 
However our interest is in using this in the three par- 
ticle case. The (a, a) element will contribute thrice to 
the (1,1) component of the three body Q-box and once 
to the (2, 2) case. Note that for the latter the use of 
— 1^ for the unperturbed energies requires that the en- 
ergy denominator be modified to 3^. This can be viewed 
as arising from diagrams with a one-body insertions of 
— ^ folded out to all orders to correct the energy of the 
model space -I- state. The two-body part of the (5, b) 
element contributes twice to the (2,2) component; the 
one-body part is associated with the single -I- state and 
therefore only contributes once. Finally the {b, a) element 
contributes thrice to the (2, 1) component multiplied by 

a factor due to the normalization in Eqs. H10|l and 
()24|l . In this way the one- plus two-body contribution to 
the three-body Q-box is found to be 

31/2 

(i|g|i) = -— , 



(2|Q|1) = (1|Q|2) = V3J7- 
C/2 



V3VU 



2e 



(2|Q|2) = ^ + 2W- 



2e 



3C 



(26) 



Note a subtlety here. A diagram of the type shown in 
Fig. Elmight appear to be a two-body contribution to the 
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0.5 I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




FIG. 7: Eigenvalues (in units of 5) obtained with the two- 
body eflective interaction (dashed curves) compared with the 
exact values (solid curves), as a function of w/^. 

three-particle interaction. However if the non-interacting 
line is erased the remainder cannot contribute to the two- 
body interaction because the intermediate state is in the 
model space. This is no longer the case with three par- 
ticles present where the intermediate state lies outside 
the model space, thus it is properly counted as a three- 
body interaction. Comparing Eq. I|26|l to Eq. (|22(l we 
see that here the three-body interactions begin at second 
order due to the presence of state |2). Using the two- 
body approximation of Eq. (|26|) for the Q-box in Eq. 
(O yields the eigenvalues denoted by dashed lines in Fig. 
El Comparison with the exact result shows that the two- 
body approximation converges more rapidly in terms of 
the number of folds, but the converged result has a lower 
(upper) eigenvalue which is in error by 9% (4%); this is 
not insignificant. The point is also made in Fig. [7| where 
the same values of V/S, and C//^ are used, but the param- 
eter W/^ is allowed to vary. Again the solid line employs 
an effective interaction with one-, two- and three-body 
terms, while the dashed line is obtained when the three- 



body contributions are omitted. Clearly the accuracy of 
the two-body approximation is dependent on the precise 
values of the parameters used. In this example it can 
be quite accurate for one of the eigenvalues for particu- 
lar values of W/^, but then it is inaccurate for the other 
eigenvalue. 

IV. CONCLUDING REMARKS 

We have shown how the effective three-body interac- 
tion may be isolated from the nonfolded and folded dia- 
grams generated from a Q-box which contains one-, two- 
and three-body terms. In order to gain some insight this 
was applied to the Lipkin model where exact results can 
easily be obtained. Even for this simple example subtle 
points arose. Numerically we showed that for a system of 
three particles the three-body components of the effective 
interaction cannot be assumed to be negligible. Interest- 
ingly in our effective 2x2 matrix example the two-body 
interaction alone sometimes gave quite an accurate result 
for one of the eigenvalues, but then the other was rather 
inaccurate. 

This suggests that the role of three-body interactions 
in realistic situations deserves further study, in particu- 
lar in medium-heavy and heavy nuclei. As mentioned in 
the Introduction we hope to carry out such a study for 
the entire range of Sn isotopes using the formalism dis- 
cussed here. The difhculty in explaining the trend of the 
binding energies is a particular motivation, but more im- 
portant is the fact that three-body effects scale with the 
number of valence particles cubed, while two-body con- 
tributions scale with the square. Thus an examination of 
some thirty Sn isotopes should allow a rather definitive 
assessment of the role of three-body interactions to be 
made. 



Acknowledgments 

This work was supported in part by the US Depart- 
ment of Energy under grant DE-FG02-87ER40328 and 
the Research Council of Norway. 



[1] A. Polls, H. Miither, A. Faessler, T. T. S. Kuo, and E. 

Osnes, Nucl. Phys. A401, 124 (1983). 
[2] H. Miither, A. Polls, and T. T. S. Kuo, 

Nucl. Phys. A435, 548 (1985). 
[3] P. K. Rath, A. Faessler, H. Miither, and A. Watt, 

Nucl. Phys. A492, 127 (1989). 
[4] D. J. Dean, T. Engeland, M. Hjorth- Jensen, M. P. Kar- 

tamyshev, and E. Osnes, Prog. Part. Nucl. Phys. 53, 419 

(2004). 

[5] A. P. Zuker, Phys. Rev. Lett. 90, 042502 (2003). 
[6] M. Homna, T. Otsuka, B. A. Brown, and T. Mizusaki, 
Phys. Rev. C 69, 034335 (2004). 



[7] T. Otsuka, M. Homna, T. Mizusaki, N. Shimizu, and 

Y. Utsuno, Prog. Part. Nucl. Phys. 47, 319 (2001). 
[8] B. A. Brown, Prog. Part. Nucl. Phys. 47, 517 (2001), and 

references therein. 
[9] M. Hjorth- Jensen, T. T. S. Kuo, and E. Osnes, 

Phys. Rep. 261, 125 (1995). 
[10] M. Lipoglavsek et al, Phys. Rev. C 65, 021302(R) 

(2002). 

[11] M. Lipoglavsek et al, Phys. Rev. C 65, 051307(R) 
(2002). 

[12] M. Lipoglavsek et al, Phys. Rev. C 66, 011302(R) 
(2002). 



9 



[13] T. Engeland, M. Hjorth- Jensen, and E. Osnes, 

Phys. Rev. C 61, 021302(R) (2000). 
[14] J. J. Ressler, W. B. Walters, C. N. Davids, D. J. Dean, 

A. Heinz, M. Hjorth- Jensen, D. Seweryniak, and J. Sher- 

gur, Phys. Rev. C 66, 024308 (2002). 
[15] A. Holt, T. Engeland, M. Hjorth- Jensen, and E. Osnes, 

Nucl. Phys. A634, 41 (1998). 
[16] S. C. Pieper, V. R. Pandharipande, R. B. Wiringa, and 

J. Carlson, Phys. Rev. C 64, 014001 (2001). 
[17] S. C. Pieper, K. Varga, and R. B. Wiringa, Phys. Rev. C 

66, 0044310 (2002). 
[18] R. B. Wiringa and S. C. Pieper, Phys. Rev. Lett. 89, 

182501 (2002). 

[19] P. Navratil and W. E. Ormand, Phys. Rev. Lett. 88, 
152502 (2002). 

[20] P. Navratil and W. E. Ormand, Phys. Rev. C 68, 034305 
(2003). 

[21] T. T. S. Kuo and E. Osnes, Folded- diagram Theory of 
the Effective Interaction m Atomic Nuclet, Springer Lec- 



ture Notes in Physics (Springer, Berlin, 1990) Vol. 364; 
T. T. S. Kuo, Lecture notes in Physics; Topics in Nuclear 
Physics, eds. T. T. S. Kuo and S. S. M. Wong, (Springer, 
Berlin, 1981) Vol. 144, p. 248. 

[22] I. Talmi, Simple Models of Complex Nuclei, Contempo- 
rary Concepts in Physics (Harwood Academic Publish- 
ers, Chur, Switzerland, 1993) Vol. 7. 

[23] K. Kowalski, D. J. Dean, M. Hjorth- Jensen, T. Papen- 
brock, and P. Piecuch, Phys. Rev. Lett. 92, 132501 
(2004). 

[24] H. J. Lipkin, N. Meshkov, and A. J. Click, Nucl. Phys. 
62, 188 (1965). 

[25] K. J. Abraham and J. P. Vary, preprint |nu'cl-th/0410037l 
[26] B. H. Brandow, Rev. Mod. Phys. 39, 771 (1967). 
[27] P. J. Ellis and E. Osnes, Rev. Mod. Phys. 49, 777 (1977). 
[28] K. Suzuki, R. Okamoto, P. J. Ellis, and T. T. S. Kuo, 
Nucl. Phys. A567, 576 (1994). 



